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If the color Coulomb potential is confining, then the Coulomb field energy of an isolated color charge is 
infinite on an infinite lattice, even if the usual UV divergence is lattice regulated. A simple criterion for Coulomb 
confinement is that the expectation value of timelike link variables vanishes in Coulomb gauge, but it is unclear 
how this criterion is related to the spectrum of the corresponding Faddeev-Popov operator, which can be used to 
formulate a quite different criterion for Coulomb confinement. The purpose of this article is to connect the two 
seemingly different Coulomb confinement criteria, and explain the geometrical basis of the connection. 
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I. INTRODUCTION 



The Coulomb potential in non-abelian gauge theories is of 
interest for several reasons. First of all, since a confining 
Coulomb potential is a necessary (though not sufficient) con- 
dition for having a confining static quark potential [1], an un- 
derstanding of the former type of potential could be helpful in 
understanding the latter, at least in the Casimir-scaling regime. 
Secondly, the Coulomb potential may be useful in various 
hadron phenomenology and spectrum calculations, perhaps 
along the lines suggested by Szczepaniak and co-workers [2]. 
Finally, the confining Coulomb potential is an important in- 
gredient in the "gluon-chain" model [3], which is a theory of 
the formation of color electric flux tubes (for a recent devel- 
opment, see [4]). 

If there is a confining Coulomb potential, then the Coulomb 
energy of an isolated color charge in an infinite volume must 
be infinite, even with a lattice regulation of the usual ultravio- 
let divergence, and this condition can be expressed in two very 
different ways. The first, derived in ref. [5], is obtained from 
the expectation value of the non-local term in the Coulomb 
gauge Hamiltonian, and is a condition on the density of near- 
zero modes of the Faddeev-Popov operator. The second is 
the criterion that the expectation value of timelike link vari- 
ables vanishes in Coulomb gauge [6]. These two criteria are 
so different that they appear to have, at best, only a very indi- 
rect relationship. In this article I will show how one criterion 
depends on the other, and discuss the relevant properties of 
gauge orbits which underlie this dependency. 

It is important to first understand where these two differ- 
ent Coulomb confinement criteria come from. Consider the 
physical state 



(1.1) 



in Coulomb gauge, where q" is a heavy quark operator, and 
^0 is the ground state. Let 



T = exp[— (// — Eo)a] 



(1.2) 



represent the lattice theory transfer matrix, divided by a factor 
of exp(— iiofl), where Eq is the vacuum energy and a is the 



lattice spacing. Then 

exp[-E,eif] 



(Tr[f/o(x,f)]) 




(1.3) 



where E^eif is the Coulomb energy, in lattice units, of the iso- 
lated heavy quark state, and the last line holds if this self- 
energy is infinite in an infinite volume (with the UV diver- 
gence controlled by the finite lattice spacing). This means that 
Coulomb confinement, on the lattice, is equivalent to having a 
vanishing expectation value for the trace of timelike link vari- 
ables. It was noted long ago that Coulomb gauge does not 
fix the gauge completely; there is a remnant gauge symmetry 
which depends on time but is homogenous in space: 



U,{x,t) ^ g{t)U,{x,t)g\t) 
Uo{x,t) -> g{t)Uo{x,t)g\t+l) 



(1.4) 



If remnant gauge symmetry is unbroken, then (Tr[?7o]) = 
[7], and the Coulomb energy of an isolated color charge is 
infinite [6]. 

On the other hand, the Coulomb energy of an isolated 
charge can alternatively be expressed in terms of the inverse 
Faddeev-Popov operator: 



{M-^{-V^)M-^)l^ (1.5) 



where Cf is the quadratic Casimir of the fundamental repre- 
sentation of the SU(N) gauge group, g is the gauge coupling, 
and 



(1.6) 



is the Faddeev-Popov (F-P) operator, with D'j^' the covariant 
derivative. Let {A-„,(|)J^(x)} denote the set of eigenvalues and 
eigenstates of the F-P operator 



M"^<(y)=X„C(j 



(1.7) 



Then it is fairly straightforward to show that in an infinite vol- 
ume, where the spectrum of M is continuous, we have [5] 



X2 



(1.8) 
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where p{X) is the density of eigenvalues, normalized to unity. 
Coulomb confinement (Eseif = °°) requires that integral in 
(1.8) diverges from a singularity in the integrand at the lower 
limit, due to the near-zero eigenmodes. The existence of such 
eigenmodes implies that the Coulomb gauge lattice configu- 
ration, in an infinite volume, lies at the Gribov horizon, in ac- 
cordance with the Gribov-Zwanziger scenario [8]. However, 
while proximity to the horizon is certainly a necessary condi- 
tion for Coulomb confinement, it is not sufficient. The density 
and behavior of near-zero eigenmodes must also be such that 
the integral in eq. (1.8) is divergent. 

We therefore have two, apparently quite distinct, criteria for 
Coulomb confinement: 

A) lun- ^>0 (1.9) 

B) (Tr[f/o(x,f)])=0 (1.10) 

These conditions look very different. They are, in fact, as- 
sociated with two different ways of computing the Coulomb 
potential. The first is to calculate the potential 



Vc{R) = -g'Cp 



1 



[M-i(-V2)M-i] 



(1.11) 



directly, on the lattice, and a number of authors have followed 
this approach [9]. The published results provide data for the 
Coulomb potential in momentum space, which appears to go 
as l/k'^ at small k. A second, and computationally much sim- 
pler method was suggested in ref. [10], and only requires com- 
puting the correlators of timeUke link variables, in Coulomb 
gauge, at equal times. Define 



(1.12) 



and let at = t,a be the lattice spacing in the time dkection. 
Then 



exp[-^£(/?) 



7t 



= (Tr[f/J(x,f)f/o(y,f)]) (1-13) 

so that the Coulomb energy E (R) associated with a static qq 
state, in units of the spatial lattice spacing a, is given by the 
timelike link-link correlator 



(Tr[{/o^(x,Ot/o(y,0]) 



= V'c(^)+ const 



(1.14) 



where R= |x — y|. This approach has been followed (at ^ = 1) 
in refs. [10, 11]. The Coulomb string tension Oc is extracted 
from the exponential falloff of the timelike link-link correlator 



(Tr[f/J(x,f)t/o(y,0]) 



(1.15) 



and of course this exponential falloff is only possible if condi- 
tion B (eq. (1.10) above) is satisfied. The interesting question 



is why (Tr[C/o]) = 0, and how this condition is related to the 
Coulomb confinement criterion A (eq. (1.9)), which is formu- 
lated in terms of the spectrum of the F-P operator. 



II. GAUGE ORBITS AND THEIR NEAR-TANGENTIAL 
INTERSECTIONS 

The suggestion 1 will make here is that the condition 
(Tr[C/o]} = 0, and the existence of a finite correlation length 
among timelike link variables, is associated with the way 
in which a typical gauge orbit intersects the submanifold of 
gauge fields satisfying the Coulomb gauge condition, and that 
this in turn is related to the density of near-zero F-P eigen- 
modes. In continuum notation, let A^(x) represent a gauge 
field at some fixed time t (spatial index k= 1,2,3), and 



F^[A] = V ■A"{x)=0 



(2.1) 



is the Coulomb gauge condition. The Fadeev-Popov (F-P) op- 
erator M is given by 



<[A] = -g 



where 



-3,Df5(x-y) 



^(x)=exp[/0«(x)r,] 



(2.2) 



(2.3) 



is a gauge transformation. 

Let C represent the hypersurface, in the space of gauge 
fields Aa(x) in D = 3 dimensions, satisfying the gauge con- 
dition condition (2.1). For a given A' G C, the correspond- 
ing {(])„} can be thought of as a set of orthonormal unit vec- 
tors which span the tangent space, at the point g{x) = 1, of 
the space of all gauge transformations. Since the gauge orbit 
0[A'] consists of the set of all configurations go A', the {(])„} 
also map to a set of directions spanning the tangent space of 
the gauge orbit, at the point A' G C . 

It was found in ref. [5] that for typical (i.e. Monte-Carlo 
generated) lattices, transformed to Coulomb gauge, there is a 
very large number of near-zero modes with A-„ -C 1, at least as 
compared to the corresponding spectrum of the free-field op- 
erator — 5"'' V^, and this number grows linearly with the lattice 
volume. Near zero-modes have a geometrical intepretation: 
these are "flat" directions on the gauge orbit at point A, which 
are nearly tangential to C- Said in another way: a great many 
directions on C run nearly parallel to the gauge orbits. 

To understand the implications of this fact, let us consider 
two points A^(x) and A^(x) + 5A^(x) with an infinitesimal 
separation 5A in the space of all three-dimensional gauge field 
configurations (a time slice of the D~4 configurations, ignor- 
ing the Ao component). Let gi and g2 be the gauge transfor- 
mations which bring A''^ =A andA^^' =A + 5A, respectively, 
onto C (Fig. 1) . Of course, gi and g2 are not unique because 
of Gribov copies, and also because of the remnant global 
gauge symmetry (1.4) allowed by Coulomb gauge. However, 
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gauge orbits 



where 




FIG. 1 : If the gauge-fixing hypersurface C defined by V • A = is 
almost tangential to many directions on a typical gauge orbit, then 
the gauge transformations and g2 which take two nearby configu- 
rations A''' and into Coulomb gauge may be very different. 



for a given gi the ambiguity in g2 can be eliminated by the 
requirement that g2 brings g2 {A + 5A) as close as possible 
to gioA. Then for 5A infinitesimal, the deviation of g2 from 
gi must also be infinitesimal, and we can write 



^2(x) = exp[/5e«(x)r,]gi(x) 



(2.4) 



where 71, = ^Oa for the SU(2) group used below. The devia- 
tion 89 is determined from the condition that 



F" 



g2o{A + ?>A) 







(2.5) 



Expanding this condition to first order in a functional Taylor 
series, and taking account of F[gi oA] = 0, we have 



5^" 



8e'(y)- 



8F" 



8A'l{y)=Q (2.6) 



where repeated indices are summed, repeated coordinates (y) 
are integrated, and the functional derivatives are evaluated at 
50 = 5A = 0. Then, from eq. (2.2), 



M",^[gioA]5d'\y)^gdfix) 



(2.7) 



= 5,{lTr|s,(x)<JjsI(x)<jJSAf(x)} (2.8) 

Now let {(])„} be the eigenstates of the F-P operatorM at gi oA, 
and expand ' 

;i>3 



5r(x) = E^/^cw 



;7>3 



(2.9) 



n>3 



The restriction to « > 3 in the above summations has to do 
with remnant global gauge symmetry in Coulomb gauge. For 
the SU(2) gauge group, the F-P operator has three exact zero 
modes (^i_3 = 0) on a finite volume with periodic boundary 
conditions, corresponding to the fact that if a configuration A 
satisfies the Coulomb gauge condition, so does goA for a spa- 
tially independent (i.e. global) gauge transformation. For this 
reason, the coefficients 50 1 -3 are not determined by the condi- 
tion (2.5), and can be set to anything we choose; in particular 
they can be set to zero. Note also that 5/i-3 = 0. This is be- 
cause the three zero modes are constant in space, while 5/(x) 
is a total derivative, so the inner product of ((|)„|5/) vanishes 
for« = 1,2,3. 



Substituting the expansions (2.9) into (2.7), we find for all 
n>3 that 



80„ 



(2.10) 



This equation takes us to the crux of the matter. The 5/„ coef- 
ficients depend on 5A^'(x), which is small but otherwise arbi- 
trary, so near-zero A,,, correspond in general to very large 50„. 
If there are only a few near-zero eigenmodes, then the few 
large 50„ may not contribute very much to 50 (x). On the other 
hand, if there are a large number of near-zero eigenmodes, 
then 50 (x) may also be large, and the deviation between A^'' 
and A'^' will be greatly magnified upon gauge-fixing to C, as 
illustrated in Fig. 1. 



The quantity to consider is the mean-square value of 50(x), 



The are real-valued, so the complex-conjugation symbol in eqs. (2.9), 
(2.12), and (2.15) below is superfluous. It is retained nonetheless, to indi- 
cate that the eigenstate would be a bra vector in bra-ket notation. 
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which is given by 

t2 



501 = 



V3- 



1 

1 



d\ 5e''(x)5e"(x) 



dX 



where 



8/(X)= / d^xt^l*{x)dkTi[gMk8[Oa] 



(2.11) 



(2.12) 



and V3 is the three-volume of a time-slice. Let us consider the 

magnitude of |59| for a "typical" 5A. To derive this quan- 
tity, we need to average over the 8A with some reasonable, 
gauge-invariant probability measure. ^ The simplest is a gaus- 
sian measure 

(G[5A]),„„,, 
^ JD5A e[8A]exp ^ j^/^lx 8A^(x)5A^(x)] 
/D5A exp [-^/^/^''x 8A^'(x)8A^'(x)] 

where e is an infinitesimal constant. In this measure 



(2.13) 



(5Af(x)5A5(y)) 



gauss 



:e8„W(x-y) 



(2.14) 



Then 



§6 1 )s'ni(ss 

= g2| ^^^(iCaiC (x))Trfei(x)a,^I(x)o„] 

5Af(x)Tr[gi(y)o,gi(y)a,]5A5(y)(8^;<(y)); 



eg / dX 



X2 



(2.15) 



The above expression still depends on the initial choice of A, 
since the F-P eigenmodes are evaluated at gi oA, but if we 
now take the vacuum expectation value, then the integral is 
precisely the same as the integral that appears in the expres- 
sion for the Coulomb energy of an isolated charge, shown in 
eq. (1.8).3 

There are now two possible scenarios, depending on 
whether or not the Coulomb confinement criterion A (eq. 



^ Note that SA transforms homogeneously, i.e. 8A gSAg' under a gauge 
transformation. 

^ We make use here of the fact that (Q[gi oA]), evaluated without gauge- 
fixing, is the same as evaluated in Coulomb gauge, as first pointed 
out by Mandula and Ogilivie [12]. 



(1.9)) is satisfied. Let 

G(x) = exp[/8e"(x)r„] =^2(x)^l(x) (2.16) 

serve to quantify the deviation between gauge transformations 
^1 and g2- Making use of the global remnant symmetry (1.4), 
it is always possible to set G(x) = 1 at one particular site 
X = xq. If there is no Coulomb confinement, so the integral 
in (2.15) is finite, then 59(x) is everywhere small for suffi- 
ciently small 5A. This means that G(x) « 1 and g2(x) ~ gi (x) 
everywhere in space. On the other hand, if the Coulomb con- 
finement condition A is satisfied, a quite different scenario is 
possible. It can then happen that no matter how small the mag- 
nitude of 8A, the non-compact variable 86 (x) becomes large at 
sufficiently large |x — xq | , and 89 is a random variable because 
8A is a random variable. In this case, the most likely behavior 
is that as x varies, G(x) wanders over the entire group mani- 
fold, averaging to zero in an infinite volume. Assuming that 
G(x) averages to zero, a better measure of the ^1,2 deviation is 
provided by the correlation length among the G(x), extracted 
from the correlator 



D{R) = — ( <i^.xiTr[G(x)G^(x- 
V3 J 



^R)] 



rather than G(x) itself. In particular, if D{R) goes as 



D{R) 



-fjR 



(2.17) 



(2.18) 



at large R, then there is a finite correlation length Ig = /j^ , 
and (G) =0. As 8A ^ (so g2 ~^ gi), we would expect this 
correlation length to go to infinity, i.e. /j ^ 0. Confinement 
criterion A is a necessary condition for this kind of behavior, 
but since eq. (2.6) is not necessarily valid for 89 ~ 0(1), we 
must resort, in this case, to a numerical investigation. 

To summarize: Near-zero modes of the F-P operator corre- 
spond to directions in the gauge orbit 0[A] which are nearly 
tangential to the gauge-fixed hypersurface C, at the point 
where 0[A] intersects C. Large numbers of tangential di- 
rections should have the following consequence: if A is a 
typical gauge field on a time slice, and A + 8A is a nearby 
configuration on the same timeslice, then the gauge transfor- 
mations gi and g2 which take A and A + 8A into Coulomb 
gauge will be wildly different, even for relatively small 8A. 
More concretely, in theories where the Coulomb energy of 
an isolated color charge is infinite, eq. (2.15) suggests that 
G(x) = g2{x)gl{x) is a random variable ((G(x)) = 0) with a 
finite GG correlation length. On the lattice it is possible to test 
this possibility, by calculating /)(/?) numerically. 



III. NUMERICAL RESULTS 

We begin by generating, via lattice Monte Carlo simula- 
tions, a set of thermalized SU(2) lattices, using the usual Wil- 
son action in D = 4 dimensions at p = 2.2. Take any time 
slice of such a lattice, and denote the link variables as uj,^\x). 
The configuration is fixed to Coulomb gauge by the over- 
relaxation method, and the gauge transformation (a product 
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of the transformations obtained at each over-relaxation sweep) 
taking C/'^^ to Coulomb gauge is denoted (x). Next we con- 
struct a "nearby" lattice U^^\x) by adding a small amount of 
noise to each link variable in the original (non-gauge fixed) 
t/(i),i.e. 



Ul'\x)=p,{x)ur{x) 



(1) 



(3.1) 



where pir(x) is a stochastic SU(2)-valued "noise" field, biased 
towards the identity, and generated independently at each link 
with probability distribution 



prob. measure o= exp 



KiTr(p) 



dp 



(3.2) 



where dp is the Haar measure. With this probability measure, 
the average value of Tr(p) as a function of K is 



Tr[p] 



/i(K) 



(3.3) 



Having generated C/*^^' in this way, we fix it to Coulomb 
gauge by the same over-relaxation procedure that was ap- 
plied to and obtain g2(x). From this we construct 
G(x) = g2{x)g[ (x), and obtain, on a lattice of extension L in 
the spatial directions, the correlator 



D{R) = jjY,Ti[G{x)G'{x + R)] 



(3.4) 



The final step is to average the values of D{R) obtained on 
every time-slice of every lattice of a set of independent, ther- 
malized lattices. The result, for average Tr(p)=0.75 (k = 5.67) 
and a variety of L"* lattice volumes, is shown in Fig. 2. It is 
clear that D{R) does indeed fall off exponentially, with inverse 
correlation length /j = 0.30. 

As K — > oo and Tr(p) ^ 1, it must be that the inverse cor- 
relation length jj goes to zero. Fig. 3 is a plot of /j vs. aver- 
age Tr(p). For lattice configurations U^^^ generated (without 
gauge fixing) at p = 2.2, the results are consistent with 



A< = c(l-iTr[p] 



(3.5) 



and c = 1.23, as the average Tr[p] 1. On a finite lattice (the 
maximum size used here is = 22''^), the practical constraint 
on K is that it should not be so large that finite size effects 
are dominant. Note also that g2 gi only for C/'^' approach- 
ing f/''), rather than approaching some arbitrary gauge copy 
of C/f''. If, e.g., f/'^' is taken to be a random gauge copy 
ofU^^K and the f/C'^) 

are transformed via over-relaxation to 
Coulomb gauge, then the gauge-fixed configurations are, in 
general, Gribov copies of one another. Numerically it is found 
that D{R) is consistent with zero in this case, for all R> 0. 

Finally, since the conjecture is that (G) = results from a 
high density of near-zero F-P eigenvalues, as compared to the 
density in a free field or non-confining theory, we should find, 
conversely, that (i) (G) ^ 0; and (ii) D{R) has a non-zero limit 
as /? oo; when G(x) is evaluated for field configurations [/* ' ' 



<G G> Correlator, pure gauge theory 
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FIG. 2: The gauge transformation correlator D{R) vs. R. The cor- 
relator is evaluated, at various L"* lattice volumes, for configurations 
{/''^ generated by lattice Monte Carlo in SU(2) pure gauge theory at 
P = 2.2, and C''^^ derived from c''' with noise parameter K = 5.67 
(^Tr(p) = 0.75). The straight line is a fit of exp{~^R) to the data at 
low R. 
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1 - Tr(p)/2 

FIG. 3: Inverse correlation length /j, extracted from D{R) calculated 
for various values of the noise parameter K, corresponding to ^Tr(p) 
in the range 0.55 — 0.85. Errorbars are smaller than the symbol size. 
The straight line is a fit to the first four data points. 



which do not have this high density of near-zero F-P eigenval- 
ues. More precisely, we expect (G) ^ for configurations in 
which the rhs of (2.11) is finite in the infinite volume limit. 
Such configurations can be generated, e.g., by center vortex 
removal in confining lattices, or alternatively by Monte Carlo 
simulation of an SU(2) gauge-Higgs theory 



5 = p ^Tr[U{P)] +y^^Tr[(pt(x)f/^(x)(p(x + ?)] 

plaq x.fi 



(3.6) 



in the "Higgs-like" region of parameter space. Here U (P) de- 
notes one-plaquette loops, and (p(x) is an SU(2) matrix-valued 
Higgs field. It was shown in [5] that the eigenvalue densities 
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<G G> Correlator, gauge-Higgs theory 



which implies Coulomb confinement? 
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FIG. 4: Same as Fig. 2, except that the t/''' are generated from lattice 
Monte Carlo simulation of SU(2) gauge-Higgs theory (eq. (3.6)), in 
the Higgs-like region at couplings P = 2.2, Y= 1-2. Noise parameter 
K = 5.67 is the same as in Fig. 2. 



are qualitatively very similar in vortex-removed and gauge- 
Higgs configurations, and are essentially a perturbation of the 
free field result. The result for D{R) in a gauge-Higgs theory 
at p = 2.2, g = 1 .2 is shown in Fig. 4. Here we see that D{R) 
does not have an exponential falloff; in fact it appears to ap- 
proach a limiting value. Hence /J = and (G) ^ in this case, 
as expected. 



IV. COULOMB CONFINEMENT 

It is time to return to the question posed in the Introduc- 
tion: how do near-zero F-P eigenmodes enforce (Tr[f/o]) = 0, 



Consider a Monte Carlo simulation carried out in temporal 
gauge, with asymmetry parameter t,. The SU(2) Wilson action 
is given by 

S = ^i:E3Tr[C/,(x,0f//(x,f + l)] 

+ P^L L iTr[t/,-(x, t)Uj{x + ei,t)uj{x + ej,t)u] (x, f )] 

(4.1) 

Suppose that ^ <C 1 . In that case, the first term in the action 
will require that link variables t4(x,f) and L'<;(x,f + 1) are 
almost identical. This means that we may think of the sets of 
link variables at fixed times t and f + 1 as being an instance of 
"nearby" D ~ 3 dimensional lattice configurations f/''^ and 
t/'^'. In fact, writing 

f/i(x,f + l)=ri^(x,f)t/i(x,f) (4.2) 
the lattice action has the form 

^ = IL L 5Ti-[m-(x,r)] + P^i:i: ^Tr[t/(P„-(x,r))] 

^ x.t k=l x,r ,</ 

(4.3) 

where U{Pij{x,t)) is a loop around the space-like plaquette 
Pij{x,t). Then, if Q is any functional of the r[k{x,t) at fixed t, 



J 



(e) = f DU{t + l)DU{t)%\u{t + l)\Q['q{t)]exp ^ iTr[Tii.(x,f)] x space-hkeplaquetti 

i X k=l 



Uit) 



I Dy\{t)DU{t) n[Tl(Ot^(0]^o[f/(f)]e[Tl(0]exp 



R 3 

— ^ ^ 5Tr[ri<-(x,f)] X space-like plaquettes 

^ X k=l 



(4.4) 



where "Polf^] is the ground state (i.e. lowest energy eigenstate 
of the transfer matrix) in temporal gauge. For ^ ^ 1 we have, 
to leading order in ^, 



(2) = j e[Ti(0]exp Iy^Y. 2Tr[il.(x, 

•' S X i:=l 



(4.5) 



So to leading order there are no correlations between the 



y\k{x,t) at different lattice sites,^ and the probability measure 
for r|i.(x) at a given link A:,x is simply 



prob. measure °^ exp 



(4.6) 



For correlations among electric field operators at different sites, and also to 
compute the Hamiltonian operator, one must of course keep tlie subleading 
terms. 
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where dr[ is the Haar measure. Identifying K = P/^, this is 
exactly the same as the probability measure for the noise vari- 
able p of the previous section, and for ^ ^ 1 it is strongly 
biased towards the identity. Therefore, f/''' = {Uk{x,t)} and 
— {Uk{x,t + 1)} are close together in the space of D = 3 
lattice configurations. 

Now to compute the Coulomb potential, we transform the 
temporal gauge lattice configurations to Coulomb gauge via, 
e.g., over-relaxation. These transformations can be computed 
independently at each time slice. Pick any time f , and denote 

f/f(x,f) = gi{x)Ukix,t)g\{x + tk) 
C/f(x,f + l) = g2ix)Ukix,t + l)glix + tk) (4.7) 

where the "C" superscript indicates the link variables trans- 
formed to Coulomb gauge, with and g2 the gauge trans- 
formations which take the temporal gauge lattices, at times t 
and t + l respectively, into Coulomb gauge. Then the timelike 
link variables at time f are simply a product of these tranfor- 
mations, i.e. 



U^ix,t) 



G\x) 



But this means that 



(4.8) 



(4.9) 



(Tr[C/o(x,0]) = (Tr[G(x)]) 
and also, referring back to eq. ( 1 . 14), that 

Vc{R) = -^"'log[(Tr[G(x)G^(x + R)])j -const. (4.10) 

The conclusion is that the property (G) = 0, which requires 
that Coulomb confinement condition A (eq. (1.9)) is satisfied, 
in turn impUes Coulomb confinement condition B, i.e. (Uq) ~ 
0. The latter condition also tells us that remnant global gauge 
symmetry in Coulomb gauge, shown in eq. (1.4), is unbroken 
(cf. the discussion in ref. [6] on this point). Moreover, a finite 
correlation length among the G(x) implies a linear Coulomb 
potential. 

In the notation of the previous section, G(x) = exp[i59(x)], 
and from eq. (4.8), we see that when G(x) is obtained from 
temporal gauge, 59" plays the role of §^q(x). It must be 
stressed, however, that while we can obtain Uq by exponen- 
tiating 59 ~ g^Ao, we cannot obtain 59 by simply taking the 
logarithm of Uq. There is an issue of which branch of the 
logarithm to choose, and choosing, e.g., the principal value 
prescription, so that 59" and run over a finite range, would 
mean that the 00 component of the gluon propagator is strictly 
bounded. In fact, in contrast to the UqUq correlator, the lat- 
tice (AqAq) correlator cannot possibly result in a confining 
Coulomb potential, as seen explicitly in the Appendix. 

For calculation of the Coulomb potential, we are consid- 
ering iy*^'^ and C/'^^ as adjacent time-slices of a thermalized 
configuration in temporal gauge, whereas in the simulations of 
the previous section, t/''^ is a time-slice of a thermalized lat- 
tice generated without any gauge fixing, and U^^^ is obtained 



from C/('^ by adding a little noise. According to the geomet- 
ric picture advocated above, this difference is unimportant so 
far as the finite correlation length is concerned. The crucial 
property leading to a finite correlation length among the G(x) 
is that f/'^^ has a high density of near-zero F-P modes when 
transformed to Coulomb gauge, and that C/'^^ is obtained by 
a small displacement from [/''^ in field space, in a random 
direction. 

Two last comments are in order, regarding the continuous 
time ^ limit. First, from eqs. (2.17) and (1.15), we see 
that in units of the lattice spacing in the spatial directions. 



(4.11) 



so if Oc is finite and non-zero in the continuous time limit, it 
must be that the inverse correlation length /j is proportional to 
^ as ^ ^ 0. Now, from the probability measure (4.6) for r|, we 
have, for small ^/p. 



;^Tr[ii]) 



or 



/2(P/^) 



l-(5Tr[Ti]) 



(4.12) 



(4.13) 



Therefore, if 0= ^ as ^ ^ 0, it must also be true that /j oc 
1 — (^Tr[r|]) in the same limit. This seems consistent with the 
data in Fig. 3 of the previous section. 

The second comment is that, since the timelike link 
correlation length Z = /j^' runs to infinity in the continuous 
time limit, there will naturally be long range correlations in 
that limit between two timelike link variables, representing a 
static quark-antiquark pair, and a timelike plaquette variable. 
What this means is that the Coulomb electric field, although 
confining, is not collimated into a flux tube. This is the 
"dipole problem" associated with all models of confinement 
based on one-particle exchange forces. In general such 
models are prone to long-range dipole fields associated with 
static charges, and long-range van der Waals forces among 
hadrons. This problem may be solved, or at least alleviated, 
in the framework of the gluon chain model (cf. ref. [4]). 



V. CONCLUSIONS 

We have seen that disorder ((Tr[t/o]) = 0) in the Coulomb 
gauge timelike Unk variables, which implies Coulomb con- 
finement, can be traced to the fact that nearby gauge con- 
figurations, when transformed to Coulomb gauge, wind up 
far apart in the space of lattice configurations. This feature 
of gauge orbits requires that a certain condition on near-zero 
modes of the Fadeev-Popov operator is satisfied, and this is in 
fact the same condition derived from requiring that Coulomb 
self-energy of an isolated color charge is infinite. 
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Faddeev-Popov near-zero modes can also be thought of as 
near-invariances of the gauge-fixing condition. The occur- 
rence of a great number of such near-zero modes means that, 
while the Coulomb gauge (in, say, the fundamental modular 
region) may indeed be a complete gauge-fixing condition, for 
typical lattices in a confining theory it is just barely so; i.e. 
there are many directions in the gauge orbit which lift the 
gauge-fixed lattice only slightly away from the gauge-fixing 
hypersurface; the gauge orbits are almost tangential, in many 
directions, to the gauge-fixing hypersurface. It is not entirely 
clear why typical gauge orbits in a pure gauge theory have this 
property, while typical gauge orbits in a gauge-Higgs theory 
do not. Proximity of a gauge orbit to the Gribov horizon is no 
doubt necessary but it is not sufficient, since it is known, e.g., 
that any lattice configuration in Coulomb gauge lying entirely 
in an abelian or center subgroup of the gauge group is on the 
Gribov horizon [6], but not all abeUan and center configura- 
tions are Coulomb confining. While there have been some 
very interesting studies relevant to the F-P eigenmode spec- 
trum (see, in particular, refs. [14] and [15]), and it is known 
that center vortex removal has a drastic effect [5], it is prob- 
ably fair to say that the density of Faddeev-Popov near-zero 
modes found in confining and non-confining theories is not 
yet well understood. 

Acknowledgments 



it may be expected that 

gooiR) = (^Tr[Ao(x,f)Ao(y,0]) 

Or 

~ -^R (A.2) 

asymptotically. This relation may be true, but the problem 
is that there is no way to verify it on the lattice. The reason 
is simple: Ao(x) is basically the logarithm of t/o(x), but the 
logarithm of Uo{x) is not unique. It is necessary to choose 
a branch of the logarithm, and we have no way of knowing 
which is the correct branch to choose. Choosing one particu- 
lar branch (e.g. via a principle value prescription) cannot pos- 
sibly result in a propagator satisfying (A.2), because R is only 
limited by the lattice size, while the magnitude of Aq, and like- 
wise Ao(x)Ao(y), is strictly bounded. 

The difficulty is quite clearly illustrated by an explicit 
Monte Carlo calculation of the lattice gluon propagator 
(}Oo{R) carried out at p = 2.2, ^ = 1 on a 16^ lattice. To 
eliminate the multi-valuedness ambiguity, we extract Aq from 
Uq = exp[;Ao • o] with the condition that < |A()| < K. The 
result is shown in Fig. 5. In contrast to the potential extracted 
from timelike link variables via eq. (1.14) [10, 11], there is no 
hint of a confining potential in the data for the gluon propaga- 
tor ^00 (^)- 
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APPENDIX: THE AMBIGUOUS AqAo PROPAGATOR 

According to eq. (1.14) above, the Coulomb field energy of 
a quark-antiquark state can be extracted from the logarithm 
of the timelike link-link correlator, and this correlator can be 
formally expressed in terms of the Aq component of the gluon 
field 



(Tr[f/J(x,OC/o(y,f)]) 



Tr 



(A.l) 



In the spirit of exponentiating ladder diagrams, it may be ar- 
gued that the exponential falloff of the link-link correlator is 
due to a confining gluon propagator, since the instantaneous 
part of the 00 component of the gluon propagator is thought 
to be proportional to the Coulomb potential [13] . Therefore, 
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FIG. 5: Monte Carlo data for the 00 component of the equal-times 
gluon propagator, with extracted from the timelike link variable as 
described in the text. There is no evidence of a confining Coulomb 
potential. The simulation is carried out on a 16'^ lattice in SU(2) 
lattice gauge theory at (3 = 2.2. 
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